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Abstract. We formulate general plant-herbivore interaction models with monotone 
plant growth functions (rates). We study the impact of monotone plant growth func- 
tions in general plant-herbivore models on their dynamics. Our study shows that all 
monotone plant growth models generate a unique interior equilibrium and they are uni- 
form persistent under certain range of parameters values. However, if the attacking rate 
of herbivore is too small or the quantity of plant is not enough, then herbivore goes ex- 
tinct. Moreover, these models lead to noise sensitive bursting which can be identified as 
a dynamical mechanism for almost periodic outbreaks of the herbivore infestation. Mon- 
tone and non-monotone plant growth models are contrasted with respect to bistability 
and crises of chaotic attractors. 

Keywords: Monotone growth models; uniformly persistent; Neimark-Sacker bifurcation; 
heteroclinic bifurcation; periodic infestations; bistability; noise bursting; crisis of chaos. 

1. Introduction 

Interactions between plants and herbivores have been studied by ecologists for many 
decades. One focus of research is the effects of herbivores on plant dynamics ^3J. In 
contrast, there is strong ecological evidence indicating that the population dynamics 
of plants has an important effect on the plant-herbivore interactions. In this article, 
we investigate how plants with different population dynamics contribute to the 
interactions. Models for plant growth vary strongly [1]: Table [T] lists eight discrete- 
time models of plant population growth. The first seven models are introduced in 
the paper by Law and Watkinson |32| without inter-specific competition. All models 
are seasonal (discrete) models of the form Pt+i — Ptf{Pt), where Pf is the density 
of a plant in season t and /(.) the per capita growth rate. In the absence of intra- 
specific competition the latter is given by /(O), i.e. I + q in models 1-3 and w in 
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models 4-8. The equilibrium density of the plant is given by K . The parameter c is 
the space per plant at which interference with neighbors becomes appreciable |32) . 
The interpretation of the power parameter, 6, depends on the model. Generally, 
these models fall into two classes, depending on whether Pf{P) is a monotone 
function of P or not. Models 1-3 are unimodal, i.e., they have a single hump. They 
lead to complicated dynamics including period doubling, period windows and chaos 
[8]. Models 4-8 are monotone, leading to much simpler dynamics. Model 8 has a 
growth function of Holling-Type III [22j . 



Table 1. Growth models of plant population density 
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Notice that Model #2 is the well known Ricker model [53] which is unimodal and 
usually written as 

Pt+i ^Pte'-^i-*) (1.1) 
while Model #4 is the Beverton-Holt model [2] usually written as 

e ^ K + Pt{l - e '^j 



The dynamics of the Ricker model (1.1) has been well studied. It shows period- 
doubling, chaos and period windows. A plant-herbivore model with Ricer dynamics 
in plant has been studied in [15] (also see similar models in [TT] and [19]) showing 
many forms of complex dynamics. 

There are fair amount of literatures on seasonal (discrete) multi-species interac- 
tion or stage structure models (e.g., Abbott and Dwyer [T7]; Gushing [5]; Dhirasak- 
danon [7] ; Jang [2] ; Kang et al [TS] ; Kon [T^] ; Roeger [5S] ; Salceanu [33] ; Tuda and 
Iwasa [30] )i among which a few studies are related to discrete prey-predator (or 
host-parasite) interaction models (e.g., Jang [HI; Kang et al [TS]; Kon [IH])- Tuda 
and Iwasa [3D] developed scramble-type and contest-type models to examine an evo- 
lutionary shift in the mode of competition among the bean weevils. Jang |14| studied 



January 12, 2013 12:27 WSPC/INSTRUCTION FILE monotonefinal 



Noise and seasonal effects on the dynamics of plant-herbivore models with monotonic plant growth functions 3 

a discrete-time Beverton-Holt stock recruitment model with Alice effects. Kang et 
al [T3] and Kon (THj studied a discrete plant-herbivore (or host-parasite) interaction 
model with Ricker dynamics as the growth function of plant (or host in Kon [Hj). 
In this article, we investigate the impact of general monotone plant growth models 
on the dynamics of plant-herbivore interaction. Our study is different from others 
and our results are new. We show that all monotone plant growth models generate 



a unique interior equilibrium (Theorem 4.3 ) and they are uniformly persistent (see 



related definitions in [27j) for certain range of parameters values (Theorem 4.5). If 



the attacking rate of herbivore is too small or the quantity of plant is not enough. 



then herbivore goes extinct (Theorem 4.2 1. In addition, our numerical simulations 



suggest that these models lead to noise sensitive bursting which can be identified as 
a dynamical mechanism for almost periodic outbreaks of the herbivore infestation. 

The rest of paper is organized as follows. In Section 2 we define two classes 
of monotone dynamics of single plant species. In Section 3 we formulate general 
plant-herbivore models for the plant dynamics introduced in Section 2. In Section 
4 we analyze the dynamic behavior of these two general models, e.g. the global 
stability of the boundary equilibrium and uniform persistence of these models. In 
Section 5 we apply the theoretical results from Section 4 to a Beverton-Holt model 
and a HoUing-Type III model. The analysis and numerical simulations suggest that 
Beverton-Holt model goes through Neimark-Sacker bifurcation with unique periodic 
orbit for a certain set of parameters values; while Holling-Type HI model goes 
through heteroclinic bifurcation for a certain set of parameters values. Our study 
also shows that noise is an important factor for outbreak of herbivore. Finally, we 
compare monotone plant growth models to unimodal and multimodal plant growth 
models regarding their influence of plant-herbivore dynamics. 



2. Monotone growth dynamics for a single plant species 
Consider 

Pt+i = Ptfir, Pt) = F(r, Pt), t > 0. (2.1) 

where Pt is the density of biomass in plant at generation t; F{r,Pt) is the growth 
function of biomass density and /(r, Pt) is the per capita growth rate of the biomass 
density. Without intra-specific competition, we have /(r, 0) = r, i.e., r is the maxi- 



mal per capita growth rate of the plant. This simple formulation (2.1 1 can give rise 
to a great diversity of dynamical behavior, depending on the expression used for 
the growth function /(r, •) and the values given to the parameters of that function. 
Several different functions have been considered. See |S] for a partial list of models 
with per capita growth rates that decline with increasing population density: 

^^^'l^ < 0, P > 0. (2.2) 

In biological terms, this means that the per capita growth rate /(r, P) decreases due 
to negative density dependent mechanism such as intra-specific competition between 
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individuals within a population. For convenience, we use F{P),f{P) instead of 
F{r, P), f{r, P) since r is a fixed parameter. The well known prototypes of the 
model ( |2.1[ ) under this biological assumption are the Beverton-Holt and Ricker 
models. The dynamics of Ricker model has been extensively studied (e.g., [23]; (15) : 
[in]). Here, we focus on the Beverton-Holt prototype, i.e., the dynamics of the plant 
is monotonically increasing, 

F\P) = ^L^li^ > 0, P > 0. (2.3) 
dP 

We can characterize the growth models of a single plant with assumption HI or 
H2 or both HI and H2: 

HI: F(0) = 0, F(P)|p>o > 0, F'{P) > 0, and limp^+oo F{P) = 0. 
H2: /(P)|p>o > 0, /'(P) < 0, and limp^+oo /(P) = 0. 

In the biological sense, HI implies that the population density in one year is a in- 
creasing function P(P) of the density in the previous year and its per capita growth 
function /(P) may be increasing or decreasing or both, which implies that plant 
suffers from the extremes of contest intraspecific competitive interaction (Henson 
and Gushing [10 ); H2 implies that the per capita growth function of the plant is 
a decreasing function due to intra-specific competition and the population density 
of a plant can be an increasing or decreasing function or both with respect to its 
density, which implies that plant suffers from the extremes of scramble intraspecific 
competitive interactions (Henson and Gushing fTO'). In this article, we study the 
population dynamics associated with plants that satisfy HI or H2 or both HI 
and H2. The specific assumptions will be addressed in the models. The following 
proposition summarizes the dynamics of plant in the absence of herbivore. 

Proposition 2.1. 

(1) Assume that HI holds and there are n + 1 consecutive, distinct and non- 
degenerate solutions P*, i — 0,1, n of P — F(P) with the following property 

= P" < pi < • • • < P" . 

// P° is stable (unstable) then the even P* are stable (unstable) while the odd 
P* are unstable (stable). In particular, P" is always stable. Moreover, define 
the map Pt+i — F{Pt), then for any e > 0, there exists N large enough, such 
that for all t > N , we have 

Pt+i ^ F{Pt) < P'^ + e. 

(2) Assume that H2 holds, then P — Pf(P) has at most two roots, i.e., P — and 
the possible root of 1 = f{P). 

Proof. : Possible configurations of the staircase diagrams of Fig. [T] shows the al- 
ternating stable and unstable equilibria. If HI holds, then we have 

< lim F(P) = C < 1 and < P'(P") < 1. 

P— >-+oo 
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This implies that the largest equilibrium P" is locally stable. 




Fig. 1. Possible configurations of the staircase diagrams 

If an initial condition Pq satisfies Pq < P", then HI implies that 

Pi = F{Po) < F{P'') = F" 

and by induction, Pt < P" for all t > 1. In the case that the initial condition is 
larger than P", i.e., Pq > P", then HI and the fact that P" is the largest positive 
root of P = P(P) indicate that P(P) < P for ah P > P". Thus, we have follows 

pn ^ p(^pn^ < Pi = P(Po) < Pq. 

Therefore, by induction, we know that the sequence {Pt}^o decreasing and con- 
verges to P" a.s t — >■ c». This indicates that for any e > 0, there exists N large 
enough, such that for all t > N, we have Pj+i = F(Pt) < P" + e. In other cases, 
we have Pt+i = P(Pt) < P" for all t > 0. 

Since /(P) is a differentiable and strictly decreasing function of P, thus 1 = /(P) 
has at most one solution. Therefore, the statement holds. n 

3. Plant-herbivore models 

Insect and plant survival rates often appear to be non-linear functions of plant 
and insect density, respectively ([3]; [3]). In our discrete-time models, we therefore 
assume that the plant population growth is a non-linear function of herbivore and 
plant density, and that plant population growth decreases gradually with increasing 
herbivore density. Similarly, we assume that the density of herbivore population 
depends on both the plant and herbivore's density rather than only the herbivore 
density '3]. A final key feature of many plant-herbivore interactions is that, in the 
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absence of the herbivore, we have a monotone growth dynamics as discussed in the 
previous section. 

Let Pt represent the density of edible plant biomass in generation t and Ht 
represent the population density of herbivore. The effect of the herbivore on the 
plant population growth rate is described by the function ^(a, Ht) with g{a, 0) — 1. 
Here the parameter a measures the damage caused by herbivore, e.g., feeding rate. 
We assume that the herbivore population density is proportional to a function of 
plant density h{Pt) and a non-linear function of herbivore density l{Ht). Therefore, 
the structure of our models is: 

Pt+i^Ptf{Pt)9{a,Ht) (3.1) 
Ht+i = h{Pt)l{Ht) (3.2) 

Many consumer-resource models assume a non-linear relationship between resource 
population size and attack rate ([T]; [12] )■ For plants and insect herbivores, we 
similarly expect a non-linear functional relationship, due to herbivore foraging time 
and satiation. The relationship is expressed in terms of plant biomass units rather 
than population size, because herbivores are unlikely to kill entire plants ([HI; |3]). 

Our model has the following features: Without the herbivore, we assume a mono- 
tone growth rate, i.e., HI holds. The growth function F{Pt) determines the amount 
of new leaves available for consumption for the herbivore in generation t. We assume 
that the herbivores search for plants randomly. The area consumed is measured by 
the parameter a, i.e., a is a constant that correlates to the total amount of the 
biomass that an herbivore consumes. The herbivore has a one year life cycle, the 
larger a, the faster the feeding rate. After attacks by herbivores, the biomass in the 
plant population is reduced to 

Pt+i = Ptf{Pt)e-^"' (3.3) 



where g{a,Ht) in (3.1) is defined as 

g{a,Ht)^e-''"\ (3.4) 



The term h{Pt) in (3.2 1 describes how the biomass in the plants is converted to 
the biomass of the herbivore. It differs depending on the relative timing of herbivore 
feeding and growth. If the herbivore attacks the plant before the plant grows, then 
we have h{Pt) = Pt, otherwise, h{Pt) = Ptf{Pt)- Since the biomass of herbivore 
comes from whatever they eat, h{Pt) is the available biomass of a plant that can 
be converted into the herbivore's biomass. The term l{Ht) describes the fraction 
of h{Pt) that can be used by the herbivore, i.e., l{Ht) = 1 — e~°^'. Therefore, the 
evolution of the plant-herbivore system is either described by Model I : 

Pt+i = F{Pt)e-^''^ (3.5) 
Ht+i - Pt [1 - e-'^^'] (3.6) 
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describing the dynamics of a system where the plant is attacked before it has a 
chance to grow while Model II : 

Pt+i - F{Pt)e-'^"^ (3.7) 
i/t+i = F{Pt) [1 - e-'^"^] (3.8) 
describes the dynamics when the plant grows first before being attacked. 



4. Mathematical analysis 



First, we can easily see that M?, is positively invariant for both Model I and II. In 



addition, 

Lemma 4.1. If HI holds, then limsupj_j.oQ maxjPt, < P" for both Model I 
and II. 

Proof. : For Model I, 

ift+i = Pt [1 - e-''"*] < Pt 

For Model II, 

Ht+i = PtfiPt) [1 - e--^^'] < F{Pt) 



Since Condition HI holds for P(P), then from Proposition (2.1), we can conclude 
that for any e > 0, there exists N large enough, such that for all t > N, the following 
holds 

Pt+i = F{Pt)e-''"^ < F{Pt) < P" + e. 
Therefore, we have limsupt,^^^ max{Pf , Ht} < P" for both Model I and II. □ 

4.1. Equilibria and their stability 

If, in the absence of the herbivore there exist n+1 equilibria of the plant dynamics, 
then both Model I and II have n + 1 boundary equilibria of the form 



Poo = (0,0) andP.o = (P\0),z = 1,2 



, ^, . . , 1 



Their local stability can be determined by the eigenvalues of their Jacobian matrices. 

It is easy to check that the Jacobian matrices of Model I and II are identical 
at these boundary equilibria: the eigenvalues of their associated Jacobian matrix at 
(0, 0) are P'(0) and 0; the eigenvalues of their associated Jacobian matrix at (P*, 0) 
are P'(P*) and aP*. The following theorems summarize the global dynamics: 

Theorem 4.2. Assume that HI holds for both Model I and II. If F'{0) < 1 and 
(0, 0) is the only boundary equilibrium, then Model I and II are globally stable at 
(P°,0) = (0,0). More generally, if aP" < l,n E Z+, then \imt^^ Ht = for both 
Model I and II. 
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Proof. From Lemma ( 4.1 ), we know that for any e > 0, there exists N large enough, 
such that for aU t > iV, we have 

Pt+i = F{Pt)e-'"^ < F{Pt) < P" + e (4.1) 

Since aP" < 1, then for e small enough, we have 

aPt < a{P" + e) < 1 and aF{Pt) < a(P" + e) < 1 for all t > N. 

Thus, for Model I, 

Ht+i = Pt [1 - e-"^*] = HtPt^ i < aHtPt < a{P'" + e)Ht (4.2) 

tit 

and for Model II, 

i/t+i = F{Pt)Ht^ i < aHtF{Pt) < a{P^ + e)Ht (4.3) 

tit 

Therefore, we have Ht < [a(P" + ej^-^HN, for all t > N. This indicates that 
lim4_j.oo Ht = for both Model I and II. Hence solutions of Model I and II are 
globally attracted to the boundary dynamics. □ 



Theorem 4.2 indicates that herbivore can not maintain its population if its at- 
tracting rate is too small or there is no enough food, i.e., aP" < 1. In addition, the 
special case of Theorem |4.2| when n = 1 leads to the following remarks. 

Remark: Assume that the hypotheses of Theorem 4.2 hold. If n = 1, then from 



Proposition (2.1), we have 



(1) If P'(0) < 1, then Pi is a source; 

(2) If P'(0) > 1, then Pi is a sink. 

Hence, if P'(0) > 1 and n — 1, then (P^,0) attracts all nontrivial solutions. 
4.2. Unique interior equilibrium 

Interior equilibria are determined by the intersections of the nuUclines. Notice that 
if HI holds, then y = F{P) is a differentiable and monotone function of P and 
maps M+ to [0, C). Its inverse exists and can be written as P = F~^{y) which maps 
[0, C) to Similarly, if H2 holds, then y = f{P) is a differentiable and monotone 
function of P and maps M"*" to [0, M). Its inverse exists and can be written as 
P = /-i(y) which maps [0, M) to M+. Here C = P(oo) and M = /(O) are some 
positive constants. If (P, H) is an interior equilibrium, then it is the solution of the 
two equations 

(1) For Model I: 

P = r\e^") (4.4) 

P= (4.5) 

1 - e-"" ^ ' 
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(2) For Model II: 

P = (4-6) 

If F{P) is monotonically increasing, i.e., HI holds, then F^^{ ^_^_^„ ) is an increas- 
ing function of H which attains its minimum at iJ = 0, i.e., 

inl^"' { "-^\]=F-\ n = ^"'(-)- 



mm ■ 



Similarly, if f{P) is monotonically decreasing, i.e., H2 holds, then / ^(e°^) is a 
decreasing function of H which attains its maximum at iJ = 0, i.e.. 



Theorem 4.3. 

a) Assume that both HI and H2 hold for Model I, then Model I has at most one 
interior equilibrium which occurs when /~^(0) > ^. The interior equilibrium 
emerges generically through a transcritical bifurcation from the largest boundary 
equilibrium P" when P" — ^, where n > 1. 

b) Assume that HI holds for Model II, then Model II has at most one interior 
equilibrium which occurs when P^^(^) < ^- The interior equilibrium emerges 
generically through a transcritical bifurcation from the largest boundary equilib- 
rium P" when P" = - , where n > 1. 

Proof. : The proofs for a) and b) are similar. We show case b): The interior equi- 



libria of Model II are determined by the intersections of the nuUclines (4.6) and 



(4.7|. Since (4.6 1 is a decreasing function and (4.71 is an increasing function, they 



have only one interior intersection if the following inequality holds 

minjP^i I ^^-^ I I < maxj^f^^^ — | ^ P^V-) < ^. 

H>o^ \l-e-''"J' H>o^e'^^-l' a a 

The Jacobian matrix of Model II evaluated at the boundary equilibrium (P* , 0) 

is 



J\(p 



P'(P') -aP' 
aP' 



(4.8) 



with its eigenvalues as Ai = aP* and A2 = F'{P^). Thus, at the largest boundary 
equilibrium (P"',0), we have 



(1) Ai|(p,.,o) = aP" and A2|(p..o) = F'(P") < 1; 

(2) l^|(p.,o)=P"and^|(p„,o)=0 
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In the case that P" = we have Ai|(-pi 



l(P'SO) - < '^l|(P",0) - ^ "^l(P",0) 

The eigenvector associated with the eigenvalue Ai|(p,i q) = aP" is 

aP"~F'{P' 
X2 



1 and 



aP' < l,i = 1, 
9Ai, 



. . . , n — 1 and 

P" = - > 0. 
a 



-X2 



(4.9) 



If aP" = 1, then the two components of (4.9 1 have opposite signs. This implies 



that by choosing a;2 > 0, the unstable manifold of i?„o points toward the interior of 
Xii. Therefore, apply Theorem 13.5 in the book by Smoller |28| . the unique interior 
equilibrium of Model II emerges generically through a transcritical bifurcation from 
the largest boundary equilibrium P" when P" -, where n > 1. □ 



4.3. Uniform persistence of Model I and II 

We define the sets 

X = {{P,H) ■■ P >0,H > 0} 
Xn = {{P,H) eX:PH>Q} 
dXn =X\Xn 
and consider the additional hypothesis 

H3: The smallest positive root P^ of P = P(P) = P/(P) satisfies aP^ > 1, and 
in addition, /(O) > 1. 

In the following, we show that Model I and II are uniformly persistent with respect 
to {Xii, dXii) if both HI and H3 hold, i.e., for any initial condition (Pq, Hq) G Xn, 
there exists some e > such that liminft_>.oo min{P(, Ht} > e. 



Lemma 4.4. Xn anddXn are positively invariant for (3. 5)-[ 3. 6) and (3. l)-(3. 



The following theorem is the main result of this subsection. 



Theorem 4.5. IfaPi > 1, then (3.5)-{3.6) and (3.1)-(3.8) are uniformly persistent 
with respect to [Xn^dXn) provided that they satisfy both HI and 113. 



Proof. From Lemma 4.4 and Proposition 4.1 we obtain that the systems (3.5) 



(3.61 and (3.7)-(3.8l are point dissipative. This combined with the fact that the 



semiflow generated by these systems is asymptotically smooth (this is automatic, 
since the state space is in M^), gives the existence of the compact attractors of 
points for both systems (Smith and Thieme 2010 [57]). 

Notice that the omega hmit set of 5*1 = {{P,H) G : P = 0} is the trivial 
boundary equilibrium Pqo- Let i(P, H) = P be an average Lyapunov function, then 
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0. Since the system satisfies H3, then the following inequality 



we have L(P, H)\s 
holds 



Pt 



sup liminf — =sup liminf TT /(-P,)e"'''^' = sup (/(O))* > 1 



t>Q 



where {Pq, Hq) E X \ Si. Therefore, by applying Theorem 2.2 in [TT] and its corol- 



lary to the systems (3.5l-(3.6) and (3.7|-(3.8|, we obtain persistence of the plant 
population, i.e., for any initial condition Pq > 0, we have liminft_j.oo Pt ^ £■ 

The fact that the plant population is uniformly persistent implies that the system 



(3.5)-(3.6l or (3.7)-(3.8l can be restricted mXn{{P,H) G : F > e}. According 



to Proposition |2.1[ we can conclude that the omega limit sets of 52 = { {P, H) £ 
dXii : P > 0} are {EiQ, 1 < i < n}. Since aP^ > 1, Condition H3 indicates that 
aP' > 1,1 < i < n. Now define L{P,H) = as an average Lyapunov function. 



then we have L{P,H)\s2 — 0. Moreover, for the model (3.5)-(3.6l, we have 

t 



sup lim inf — - 

t>0 (Po,-H"o)-i-(P',0) Hq 



sup lim inf 

t>0 {Po:Ho)'^{P\0) 



and for the model (3.7)-(3.8l we have 



Ht 



sup liminf = sup liminf ( ]^P(Pj) 



2Hi 



t>0 (,Po,Ho)^{P\0) Ho t>Q (Po,Ho)^(P',0) 



sup (aP*)* > 1 



= sup(aP*)* > 1 
t>o 



where {Po,Hq) S Xh. Therefore, according to Theorem 2.2 and its corollary 2.3 



in [TT], we can show that the systems (3.5)-(3.6l and (3.7)-(3.8) are uniformly 
persistent. Hence, the statement holds. □ 



Remark: The arguments used to prove Theorem |4.5| are standard, which can 
be found in many literature (e.g., [T^], [T5], [TB]). 



5. Application and simulations 

5.1. The Beverton-Holt and Holling-Type III models 

In this section, we focus on two typical models for the plant dynamics and apply 
our results: 



(1) 



Pt+i = F{Pt) = 



rPt 



(5.1) 



1 + Pt 

is the Beverton-Holt model where P(Pt) satisfies the assumptions of HI and 
/(Pt) satisfies those of H2. The two equilibria are P° = and P^ = r — 1. From 



Proposition 2.1 we know that P° is a sink if r < 1; P° is a source if r > 1. In 
addition, P^, if it exists, is always a sink. 
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(2) A HoUings type III model is given by 

Pt+i = F{Pt) 



_rPl_ 
iTTf 



(5.2) 



where F{Pt) satisfies HI. The equihbria are P° = 0, = r-vj^ p2 



If r < 2 then P° is the only equilibrium and it is globally stable. If 



r > 2; P° is a sink, P^ is a source and P^ is a sink. 



The plant-herbivore models with (5.1) and (5.2) as plant dynamics become: 
Model I: 



(1) 



(2) 



Pt+i 
Ht+i 

Pt+i = 



l+Pt 



,.p2 



l+P? 



t+i 



(5.3) 
(5.4) 

(5.5) 
(5.6) 



Model II: 

(1) 



(2) 



Pt+i^ 
Ht+i — 

Pt+i = 



rPt 
l + Pt 

rPt 
l + Pt 



[1 - e-'^^*] 



rPt 

iTp? 



-aHt 



rP^ 



(5.7) 
(5.8) 

(5.9) 
(5.10) 



Applying the results of the previous section, we have the following two corollaries. 



Corollary 5.1. The three models \5. (5. 4^ , {5. 7^ - (ST^ and (5.gp-(5.iO ) have at 
most one interior equilibrium. The interior equilibria of Model I {5.3)-{5.4-) and 
Model II \5. TP - jX^ emerge through a transcritical bifurcations from the boundary 
equilibrium (P^,G) = (r — 1,0) when a{r — 1) = 1 and r > 1; The interior equilib- 



rium of Model II t57^-{5.10) emerges through a transcritical bifurcations from the 



boundary equilibrium (P , 0) = ( 



, 0) when a 



r+Vr^ — 4 



1 and r > 2. 
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Corollary 5.2. If a(r ~ 1) > 1, the systems (5.3)-(5.4-) and (5.1)-(5.8) are uni- 
formly persistent with respect to 



5.2. Periodic orbits and heteroclinic bifurcations 



The stability of the single interior equilibrium of models (5.3) to (5.10) depends on 
the values of the parameters r and a. As the values of r or a increase, the inte- 
rior equilibrium goes through a Neimark-Sacker bifurcation generating an invariant 
cycle. 

Since Model I and Model II have similar dynamics, we only focus on Model II and 
discuss the Beverton-Holt model ( 5.1 ) and the Holling-Type III model ( 5.2 ), respec- 
tively. The main difference between the Beverton-Holt model and the Holling-Type 
III model is that the Holling-Type HI model can show a heteroclinic bifurcation 
where a periodic orbit grows until it becomes a heteroclinic connection between 
boundary equilibria whereas the Beverton-Holt model does not show such a bifur- 
cation. Figure ^ shows the heteroclinic bifurcation schematically: When a = 0.71 
and r — 2.5, the system has a stable interior equilibrium (the dark dot that is in the 
middle of the figure, which is generated by the Matlab); When we increase r to 3.5 
and keep a = 0.71, the system has an invariant orbit (the grey orbit in the figure, 
which is generated by the Matlab); However, if we continue to increase the values 
of a or r, the invariant orbit disappears and the system converges to the boundary 
equilibrium (0, 0). This suggests that a heteroclinic bifurcation occurs (the dark line 
with arrows in the figure, which is generated schematically). 

Since the Beverton-Holt model only has the origin as a saddle and one other 
boundary equilibrium and since the stable manifold of the origin is the iJ-axis 
which is an invariant manifold, the periodic orbit in the interior cannot become 
heteroclinic. However, it can become very large and pass the origin arbitrarily close 
to the coordinate axes as shown in Fig. [3] Figure [3] is the numerical simulations 
generated by the Matlab for 2000 generations when a = 2 and r = 2.5, 2.7, 2.8, 3. 
When a = 2 and r = 2.5, the system has a stable interior equilibrium as shown 
in the figure (small dark dot); when r = 2,7,2.8,3, the system has an invariant 
orbit. Numerical simulations of this case hint at an interesting phenomenon: a stan- 
dard numerical simulation shows the periodic orbit disappearing and the trajectory 
approaching the nontrivial boundary equilibrium as time increases. However, that 
boundary equilibrium is a saddle and the trajectory should leave into the interior 
but it does not do so over any simulation time that we checked. The resolution of 
the puzzle comes from the considerations of the accuracy of the simulations: As 
the limit cycle gets closer to the origin the herbivore values become so small that 
they are approximated as zero. Hence the dynamics is reduced to the dynamics of 
the plant which has a stable equilibrium on the invariant manifold determined by 
H = and hence the trajectory never leaves. 

Figure |4][a) shows a "bifurcation diagram" for the Beverton-Holt model, which 
describes the Neimark-Sacker bifurcation curve (dashed line) and the "collapse 
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Plant-Herbivore model with Holling-type III model as plant dynamics 



HeterocUnic Bifurcation 
a = 0.71,r = 3.5 




Fig. 2. Schematic of tiie heteroclinic bifurcation of HoUing-Type III model happens at a = 0.71, r 
3.5. 



Plant-Herbivore model with Beverton-Holt model as plant 




0.2 0.4 0.6 0.8 1 1.2 1.4 

plant 



Fig. 3. The periodic orbit for the Beverton-Holt model when a = 2,r = 2.5, 2.7, 2.8, 3 
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curve" (solid line). The later represents an interpolation of numerical simulations 
with a and r values for which a standard matlab numerical precision simulation does 
not detect a population of the herbivore. Figure Qb) is a bifurcation diagram for 
a HoUing-Type III model showing interpolations of the Neimark-Sacker bifurcation 
curve (dashed line) and the heteroclinic bifurcation (solid line), respectively. 



I 




(a) The bifurcation diagram for the (b) The bifurcation diagram for the 
Beverton-Holt model. HoUing-Type III model. 



Fig. 4. Neimark-Sacker bifurcation and heteroclinic bifurcations 



5.3. Noise-generated outbreaks 

The extreme sensitivity of the periodic orbit in the Beverton-Holt model suggests 
that noise may play a much bigger role than previously discussed in the outbreaks 
of herbivore infestations. Once the periodic orbit disappears due to accuracy issues, 
we can make it re-appear by adding small amount of noise to the simulation: 

(1) Noise: we use positive white noise to make sure the system stays positive, i.e., 
we sample from a normal distribution but discard any negative noise sample; 

(2) Population of herbivores: for each generation, we add the noise to the herbivore, 
i.e., 

Ht+l= * „ +CjRn (5.11) 
t + -it 

where i?„ is a positive white noise as defined above and cj is the amplitude of 
the noise. See Fig. [5] for example, in this case, the amplitude of the positive 
white noise is w = 0.01. 
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At that time the trajectories look hke a randomly occurring bursting phenomenon 
that nevertheless has a well defined average periodicity (see Fig. [5]). Given by the 



Plant-Herbivore model with Beverton-Holt model as plant dynamics 

3.5 r 

a=3.95, r=4.55, w=0.01 




Fig. 5. Time series of the herbivore population for the Beverton-Holt model. The parameters are 
a = 3.95, r = 4.55 and a noise level of ui = 0.01 



exact nature of the model there will be a threshold at which the population of the 
herbivore cannot be detected in nature. We define the resident time as the time 
interval for which the population of the herbivore stays below some threshold, e.g., 
0.01 and the resident time ratio as the ratio of the residence time to the period of 
the bursting. Table [573| shows the period as a function of the mean square amplitude 
of the noise level. Figure [6] shows the resident time ratio as a function of the noise 
amplitudes. The figure is generated by calculating the resident time ratio for each 
noise amplitude for 50 trajectory with 1000 generations. The Figure shows that over 
many orders of magnitude the residence ratio stays around 80% indicating that the 
herbivore is dormant for most of the time and only appears for about 20% of its 
periodic cycle. The Table indicates that by choosing a particular noise level, we can 
control the apparent periodicity of the bursts. 



Amplitude of Noise w 


0.01 


0.001 


0.0001 


0.00001 


0.000001 


0.0000001 


Period t 


8 


10 


12 


14 


17 


19 



Table 2. Average period of the herbivore dynamics when a = 3.95, r = 4.55. 



In particular, time intervals of the herbivore outbreaks around 8—12 years can 
be generated, which fits the ecological data for gypsy moth out-breaks [H]. Also, 
for larger noise levels, the distribution of the periods is rather broad, which also 
seems to be happening for real data ([10]; [21] )■ 
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Fig. 6. The resident time ratio as a function of the noise amphtude when a = 3.95, r = 4.55 with 
a threshold of 0.01. 



6. Conclusions and additional features 

For most plant species it is conceivable that there is density dependent regulation of 
its growth. However, very few plants show periodic or strongly chaotic variation of 
the plant density from generation to generation. Hence it is important to determine 
the influence of models of monotone growth dynamics on the plant- herbivore inter- 
action model. We proved three key features of such interactions that are important 
for model building: 



• All monotone growth models generate a unique interior equilibrium. 

• Monotone growth models with just one sustainable equilibrium for the plant 
population (e.g. the Beverton-Holt model) lead to noise sensitive bursting. This 
certainly happens for many plant-herbivore systems and the dynamical mech- 
anism discussed here has not been noticed before in plant-herbivore systems 
(however, see [H]). 

• Model I and H have a uniformly persistent property if they satisfy both HI 
and H3. In particular, the Beverton-Holt model is uniformly persistent and the 
HoUing-Type III model is not. 

• The Beverton-Holt model does not have more complicated dynamics than a 
periodic orbit in the interior of the phase space. Although we cannot prove 
this, we conjecture that this is true for all models that satisfy the assumptions 
of HI and H2, i.e. have just one equilibrium for the pure plant dynamics. 



Without any claim to a complete analysis of all types of models, we note a few addi- 
tional features associated with monotone and non-monotone plant growth models. 
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(1) Bistability: The paper [TS] study plant-herbivore systems of Model II type 
with a Ricker model for the pure plant dynamics, also known as the modified 
Nicholson-Bailey model. It is shown that for a large set of parameters the sys- 
tem exhibits bistability between complicated (possibly chaotic) dynamics in the 
interior of the phase space and equally complicated dynamics on the boundary 
(Fig. [7];a)). Kon 1191 discusses a similar bistability phenomenon for Model I. 
Since unimodal maps are all topologically equivalent [S], we expect bistability 
to be a defining feature for plant-herbivore models with unimodal plant growth 
models. 

In contrast, models that satisfy the assumptions of HI and H3, e.g., the 
Beverton-Holt model cannot show bistability since the global attractor either 
is a fixed point on the boundary or some set in the interior of the phase space. 
However, it is conceivable that models that do not satisfy H3, e.g. Holling- 
Type III models, show bistability between an interior attractor and a boundary 
equilibrium that is not the largest equilibrium for the pure plant population. 

(2) Crises of Interior Attractors: All models seem to show some sort of global 
attraction to the boundary dynamics, i.e. extinction of the parasite for large 
growth rates r: 

• Unimodal models show a crisis type of bifurcation whereas the chaotic dy- 
namics in the interior collapses and the system becomes globally attracted 
to the boundary dynamics 15 . For instance the interior strange attractor 
in Fig. [7] that exists for a growth parameter of r = 3.8 will grow and hit 
the stable manifold of the boundary attractor for r = 3.85. 




Fig. 7. The interior strange attractor and the stable manifold of the boundary attractor a = 
0.95, r = 3.8. 



• HoUing-Typc III models show a heteroclinic orbit which breaks and leads 



January 12, 2013 12:27 WSPC/INSTRUCTION FILE monotonefinal 



Noise and seasonal effects on the dynamics of plant-herbivore models with monotonic plant growth functions 19 



to global attraction to a boundary fixed point. 

Conjecture: The Beverton-Holt model does not lead to complete extinction 
according to Theorem (4.5 1. However, e in Theorem (4.5) may happen 
to be small which, taking a stochastic effect into account, might lead to 
extinction of herbivores. 
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